Deformations of the Tracy-Widom distribution 
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In random matrix theory (RMT), the Tracy-Widom (TW) distribution describes 
the behavior of the largest eigenvalue. We consider here two models in which TW 
undergoes transformations. In the first one disorder is introduced in the Gaus- 
sian ensembles by superimposing an external source of randomness. A competition 
between TW and a normal (Gaussian) distribution results, depending on the spread- 
ing of the disorder. The second model consists in removing at random a fraction 
of (correlated) eigenvalues of a random matrix. The usual formalism of Fredholm 
determinants extends naturally. A continuous transition from TW to the Weilbull 
distribution, characteristc of extreme values of an uncorrelated sequence, is obtained. 

I. INTRODUCTION 

In the beginning of the 90's, Tracy and Widom (TW) derived the probability dis- 
tribution of the largest eigenvalue of random matrices belonging to the three Gaussian 
ensembles, the orthogonal (GOE), the unitary (GUE) and the symplectic (GSE)[l|. Few 
years later, Baik, Deift and Johansson proved that the longest increasing subsequence 
of a random permutation fluctuates as the largest GUE eigenvalue and triggered TW 
applications in combinatorics and other areas such as growing processes^, 0, H, |i| (see 0] 
for a review). The main ingredient in their derivation was the discovery that the formal- 
ism of random matrix theory based on Fredholm determinants and Painleve equations [s| 
which at the bulk of the spectrum is associated to integral equations with a sine-kernel, 
at the border of the spectrum it is associated to integral equations with an Airy-kernel . 
It is by now accepted that these distributions belong to universal class of extreme values 
of correlated sequences. Deviations from the TW have also been observed and studied. It 
has been found, for instance, that growing processes in random media may show regimes 
in which TW compete with other distributions [3] . 

In this general context it is important to establish links between TW and known uni- 
versal distributions of extreme values of uncorrelated sequences. In this last case, namely 
for a sequence of i.i.d. random variables, the probability that the extreme value is less 

than a given value t is exp — p(x)dx , where p(x) is their density which extends til 

x = Depending on the asymptotic behavior of the function p(x), this probabil- 

ity distribution takes three forms. With T = oo, it becomes the Gumbel distribution 
exp [— exp(— y)\ if p has a fast exponential decay and, if it decays with a power p + 1, 
it is the Frechet distribution exp(— l/y M ) (in both cases y is a properly scaled variable). 
With T finite, that is if p has a bounded support, it is the Weibull distribution exp(— y) 



where y = p(x)dx has density one. 

Largest eigenvalues of non-Gaussian ensembles (Wishart matrices) have been the sub- 
ject of recent investigations (TH] . It has been found that if the matrix elements are taken 
from a distribution with finite moments then the TW holds[ll|]. Considering instead the 
case in which the distribution of the matrix elements have long tails, it has been proven 
that when the second moment diverges, the largest eigenvalue and the largest matrix 



element follow a Frechet distribution with the same power \x < 2 [12]. This result has 



more recently been extended til \x = 4|13 



Models to describe deviations from TW have been discussed by K. Johanso nll4j| . In 
one model, he studies the behavior of the largest eigenvalue of a matrix model[15] and 
finds that it is decribed by a kernel that goes from a Poisson kernel (see Eq. f l3Tj) below) 
with an exponential density to the RMT Airy kernel. Accordingly, the distribution of 
the largest eigenvalue goes from Gumbel to TW. In another model, a deformed GUE 
ensemble is considered in which the eigenvalue density fluctuates in such a way that the 
largest eigenvalue distribution goes from TW to Gaussian. 



Following similar lines to those of [14[ the purpose of this note is to investigate other 



models that describe deformations of the TW. The first model results from superimpos- 



ing to the Gaussian fluctuations an external source of randomness |1 61] . This causes the 
eigenvalue semi-circle density to fluctuate and results in features common to growing 
processes in random media. The second model is based on the recent recognition that 
the mathematical structure of random matrix theory (RMT) also describes the statisti- 
cal properties of the eigenvalues of spectra when a fraction of eigenvalues is randomly 



removed [171 ] . As this operation reduces correlations, it describes intermediate statistics 
between RMT and Poisson statistics. Focusing on eigenvalues at the edge of the spectrum 
we show here that this leads to a transition from TW to a Weibull distribution. 
Consider the Gaussian random matrix ensembles defined by a density distribution 

faB\ f/2 

P G {H- a) = f — J exp{-a/3trH 2 ) (1) 

where f = N + /3N(N — l)/2 is the number of independent matrix elements and (3 is 
the Dyson index that takes the values 1, 2, and 4 for the Orthogonal (GOE), the Unitary 
(GUE) and the Symplectic (GSE) ensembles, respectively. In (JTj), the normalization 
constant is calculated with respect to the measure dH = Y\i dHa Ylj>i Tlk=i V^dH^. 

Let us start by recalling known facts about the eigenvalues of these Gaussian ensem- 
bles including some recent results regarding the behavior of their largest values. It is 
well known, for instance, that, to leading order, their eigenvalue density is given by the 
Wigner's semi-circle law 

rV^Na 2 - A 2 , I A |< 2aVN 



where a = 1/y/Aa is the variance of the off-diagonal matrix elements. To study the 
behavior of the largest eigenvalues in the limit of large matrix size N, one introduces the 
scaling 

A = ( 2 ^+^V (3) 
that substituted in (|2J) leads to the iV-independent density 



i v /3 s, s < 

n V ) 

0, s > 



(4) 



at the edge of the spectrum. In the scaled variable s, the probabilities E (k, s) with 
k = 0,1, 2, .. that the infinite interval (s, oo) has, respectively k eigenvalues, are obtained 
from the generating function G(s, z) through the relation [l|, [18 



G{ 8i z)=Y,{-Vr{z-VrE{n,s) 



(5) 



n=0 



such that 



E(n,t) 



dG{s,z) 
dz n 



(6) 



z=l 



For the three symmetry classes, the generating functions Gp(s,z) with /3 — 1,2 and 4 
have been derived. Starting with the unitary case, G2(s,z) can be identified with the 
Fredholm determinant associated to the integral operator acting on the interval (s, oo) 



with kernel 19 



K(x,y) 



x-y 

where Ai(s) is the Airy function. G2{s, z) is given by 

/■oo 

G2(s,z)=exp — / (x — s)q 2 (x, z)dx 
where q(s, z) satisfies the Painleve II equation 

" i o 3 

q = sq + 2q 

with boundary condition 

q(s, z) ~ V^Ai(s) when s — > oo. 
For GOE (/3 = 1) and GSE {(3 = 4) the generating functions are 20 



2 z - 1 - cosh fj,(s, z) + y/zsmhfx(s,z) 
[Gi(s, z)j = G 2 (s, z) 



and 



z-2 



[G 4 (s,z)f = G 2 (s,z) cosh 2 ^!^ 



where z = 2z — z 2 and 



q(x, z)dx. 



(7) 



(9) 



(10) 



IF 



(12) 



(13) 



We remark that the above expression for the GSE case was obtained in Refs. 19, 13 
using a scaling that assume N/2 eigenvalues. 



The above equations give a complete description of the fluctuations of the eigenvalues 
at the edge of the spectra of the Gaussian ensembles. In particular, for the largest 
eigenvalue, the TW distributions are expressed in terms of these generating functions as 

Eg,/3 (Emax < A) = Gp(s, 1). 



II. DISORDERED ENSEMBLES 

To investigate the modifications these probabilities undergo when an external source 
of randomness is superimposed to the Gaussian fluctuations, we consider disordered en- 
sembles whose matrices, if(£,a), are defined asfl6| 

where Hq is a matrix of ([I]) and £ is a positive random variable with distribution w(£) 
with average £ and variance a w . From (fl4j) and ([I]) it is deduced that the joint den- 
sity distribution of the matrix elements is a superposition of the Gaussian ensembles 
distributions weighted with u>(£), namely 

P (H; a) = J d$w{Z)P G {H; «£/£)■ (15) 

Changing variables from matrix elements to eigenvalues and eigenvectors, it is also 
found, after integrating out the eigenvectors, that the joint probability distribution of the 
eigenvalues is obtained by averaging over the joint distribution of the Gaussian ensembles. 
As a consequence, measures of this average ensemble are averages over the Gaussian 
measures. The eigenvalue density, for instance, turns out to be an average over Wigner's 
semi-circles with different radii, that is 

p(X;a) = J d£w(OV^ 2 (ON-\ 2 / [W(£)] , (16) 
where the ^-dependent variance, <r(£) is given by 

a(£) = ojlfd (17) 

The probability that the largest eigenvalue X max is smaller than a given value t can be 
calculated by evaluating the probability that the interval (t, oo) is empty. This is obtained 
by integrating the joint probability distribution of the eigenvalues in the interval (— oo, t) 
over all eigenvalues, we find 



Ep (X max < t) = J d£w(Z)E GtP [£(£, t)] (U 



with the argument of S(£, t) obtained by plugging in Eq. Q the above ^-variance, Eq. 
( JT7|) . namely 



S(£,t)=iV 1 / 6 



<K0 



(19) 



Equations (TX5T) and f[T$j) give a complete analytical description of the behavior of the 
largest eigenvalue once the function u>(£) is chosen. But even without specifying w(£), 
asymptotic results can be derived by comparing its localization, given by the ratio cr w /£, 
with that of Ec,p considered as functions of the integrand variable £ through Eq. ([19]) . 
Since the widths of these last ones depend on the matrix size N, let us introduce a positive 
parameter z such that u w /^ = N~ z which is kept fixed when the limit N — > oo is taken. 

As z > 0, when N increases, the w(£) distribution becomes more and more localized 
and, if asymptotically it can be approximated by a Gaussian, by changing the integration 
variable to 

£ = i- va w , (20) 

Eq. (TT8T) can be written as 

Ep (X max <t) = -j=j dv exp (-^j E G> p [S(v, £)] , (21) 
where the argument S(v,t), after neglecting higher order terms in 1/N, takes the form 

S{v, t) = N l ' & - vN^ 2 - z - 2^j = s - N 2 '^ z v. (22) 

In the last step of ( 1221) . Eq. (EJ was used and it was assumed that z > 2/3. Taking 
now the limit N — ► oo, the vanishing of the second term in the r.h.s. of ( 1221) makes 
S(v,t) independent of v and the distributions Ec,p(s), i.e. TW, are recovered. Values of 
the parameter z greater than 2/3 correspond to situations in which the distribution 
collapses faster than Eq q, for z < 2/3, on the other hand, it is the opposite that happens. 
To be able to get iV-independent results in this range of values of z, it is necessary to 
modify the scaling to 

s = N z ~ 1/2 (J- - 2^N^j , (23) 
in which case, the argument S \v,t(s)\ becomes 

S[v,t(s)} = N 2/3 - z (s-v). (24) 

Taking now the limit of N — » oo, E G ^{S) becomes a step function centered at v — s 
and the distribution goes to the normal distribution N(0, 1). Finally, at the critical value 
z = 2/3, from both sides, ff2Tl) converges to the convolution of the normal distribution 
and TW. In summary we have the three regimes 

E Gj/3 (s), z > 2/3 

Ep (X max <t) = \ vfe dv ex P ("t) E gA* -v), z = 2/3 (25 ) 

^ I^dv exp (-i), z< 2/3 



Consider now the case in which the distribution w(£) is independent of N. In this 
case, it is appropriate to make the parameter a equal to the matrix size N. With this 
scaling, the eigenvalues of the matrices of the average ensemble are located in the interval 
(—1, 1), and the ^-dependent argument S(£,t), Eq. ffTDl . takes the simple form S(£,t) = 



2 TV 2 / 3 {t y/^/i — 1)- This expression makes evident that when the matrix size N increases, 
for the three invariant ensembles, the function Eq^ become a step function centered at 
£ = £/t 2 . Therefore, in this regime, the probability distribution for the largest eigenvalue 
converges to 



poo 

Ep {X maX < t) = / d£w (£) (26) 

with density 

dEp (t) _ gfrgffig) 

A special choice of w(£) that already has appeared in previous studies of disordered 



ensembles [211. |22| . is that in which it is the one parameter family of Gamma distributions 

w(0 = exp(-0^-7r®. (28) 

In this case, ( I2"7j) becomes 

dEp (t) = 2^exp(-£/t 2 ) 

dt r(£)t 2f+1 ' 

that defines a long-tailed distribution of extreme values of a correlated set of points(we re- 
mark that this distribution distribution has recently been considered in random covariant 
matrices [231]). 

Notice that in f)27p the variable t is the eigenvalue itself without any edge scaling. 
This means that, in this case, fluctuations become of the order of the size of the average 
ensemble spectrum. Although (12 7p resembles a Frechet distribution, the fact that the 
power of t in the exponent is fixed at the value two makes it a different distribution. 
Nevertheless, for £ = 1 it is indeed a Frechet distribution. We remark that £ = 1 
corresponds to the critical distribution of the family defined by Eq. (128]) that separates 
the ones that converge from those that diverge at the origin. The asymptotical power- 
law decay of (1291) . similar to that of Frechet distribution, suggest that in the asymptotic 
region the extreme value behaves independently of the other eigenvalues while, in the 
internal region, the presence of the others are felt. 



III. FROM TRACY- WIDOM TO WEIBULL 

Let us now turn to the case of a model to describe largest eigenvalues of spectra in the 
intermediate regime between RMT and Poisson. This model is brought on by the fact that 
the generating function, Eq. (JS]), can be interpreted as a probability. Indeed, assuming 
with < z < 1 that the factor (1 — z) n that multiplies E(n, s) in (jSJ) is the probability 
that the n eigenvalues in the interval (s, oo) have been removed, then summing all the 
terms gives the probability that there is no level in the interval. 



A realization of this situation was considered in Ref. [24| in which the effect of re- 
moving at random a fraction 1 — / of eigenvalues of RMT spectra was investigated. In 
this case, 1 — / is the probability that a given eigenvalue has been dropped from the 
spectrum. Therefore with the identification of z with / the generating function Eq. (JSj) 



becomes the probability distribution of the largest eigenvalue for this kind of randomly 
incomplete s pec tra. 



In Refs. 

in 

|24J |. in studying the effect of incompleteness in the spectral statistics 
at the bulk, an interval of length s is increased by a factor of 1/f to compensate the 
reduction in the average number of levels inside it. Following the same idea at the edge, 
we want a scaling of the variable s such that the average number of eigenvalues in the 
interval (s, oo) remains the same when a fraction of levels is removed. This average is 

<n>=^-{-sf\ (30) 

obtained by integrating (j3J) from s to oo. Therefore, in order to keep it invariant when 
the density of eigenvalues is reduced by a factor of /, s has to be divided by f 2 ^ 3 . 

Using this scaling in the Airy- kernel, Eq. (J7J), we expect that when the limit / — > is 
taken it converges to the Poisson- kernel [3] 

(31) 

with density given by (J4j). In fact, we have 

Ai(x/ fW)M'(y/fW) - Ai(y// 2 / 3 )Ai'(s// 2/3 ) K , s 
lim — w = Kp{x, y). 32 

Denoting by e(s, f) the probability that the semi-infinite interval (s, oo) is empty, for 
the incomplete spectra it is given by 

oo 

ep(s, /) = £(!- ft E p (k, s/f 2 ' 3 ) = G p (s/f 2 ' 3 , f), (33) 

k=0 

that means that in an incomplete spectrum the largest eigenvalue can be anyone other 
of the nth largest eigenvalues. The last equality in fl33|) follows from Eq. (jSJ) and shows 
that the generating functions for the three symmetry classes contain a comprehensive 
description of the largest eigenvalues of complete and incomplete RMT spectra. 

To investigate the limit /^Owe remark that as the scaling factor f 2 / 3 appears in the 
denominator, for small / the function q(x, z) can be replaced by its asymptotic form at 
x — > ±oo. For x > 0, this is given by the exponential decay of the asymptotic behavior of 
the Airy function. So, at this positive side, q(x, z) vanishes and ip(s, f) = 1. For s < 0, 
by the same argument, the integrals can be performed from s to zero with q(x, z) replaced 
by its asymptotic expression for large negative values. For < z < 1, this expression has 
been worked out by Hastings and McLeod 125], that found 



q(x, z) ~ d(—x) sin 



^(-^) 3/2 -^ 2 log(-a;)- C 



(34) 



where d 2 = — - log(l — z). 

The integral in (jSj) contains the square of q(x, z) so, if we use the relation 2sin 2 x = 
1 + cos 2a;, we find that the generating functions have a smooth and an oscillating part. 
As the period and the amplitude of the oscillations decrease with /, in the limit of small 
/, the oscillating part averages out to zero and can be neglected. For the same reason, in 
this limit, the function fj,(s, f), Eq. ffT3"|) . vanishes. Taking then all this into account and 



after evaluating integrals, we find that when / — > 0, the largest eigenvalue distributions, 
for the three symmetry classes, converge to the Weibull distribution 



where gp = 1 for (3 = 1,2 and ^ = 1/2 for (3 = 4 (this parameter reflects the fact 
mentioned above that, in the GSE case, we are using a scaling with N/2). This Weibull 
distribution describes the extreme value of a set of uncorrelated points with a semi-circle 
density distribution. In Fig. 1, the transition from TW (/ = 1) to Weibull (/ = 0) is 
illustrated for the GUE case (J3 = 2). The oscillations at intermediate values are clearly 
seen. 



In conclusion, we have investigated the behavior of the largest eigenvalue of two models 
that generalize the RMT Gaussian ensembles. In the first one, by introducing disorder 
in the ensemble, the eigenvalue density is lead to fluctuate generating at the edge a 
competition between the Tracy- Widom distribution and the normal distribution. This 
kind of behavior has been observed in growing processes in random media. In the second 
model, the eigenvalue density is kept fixed but the correlations among the eigenvalues 
are progressively reduced. It is then observed a continuous transition from the TW to 
the Weibull distribution, characteristic of uncorrelated variables. 

This work is supported in part by the Brazilian agencies CNPq and FAPESP. 



[1] C. A. Tracy and H. Widom, Commun. Math. Phys. 159, 151 (1994) and 177, 727 (1996). 

[2] J. Baik, P. Deift and K. Johansson, J. Amer. Math. Soc. 12, 1119 (1999). 

[3] K. Johansson, Commun. Math. Phys. 209, 437 (2000). 

[4] J. Baik and E. M. Rains, J. Stat. Phys. 100, 53 (2001). 

[5] M. Prahofer and H. Spohn, Phys. Rev. Lett. 84, 4882 (2000). 

[6] J. Gravner, C. A. Tracy and H. Widom, J. Stat. Phys. 102, 1085 (2001). 

[7] C. A. Tracy and H. Widom, Proceedings of the ICM, Beijing 2002, vol. 1, 587-596. 

[8] M.L. Merita, Random Matrices (Academic Press, Boston, 1991). 

[9] J. Gravner, C. A. Tracy and H. Widom, Ann. Probab. 30, 1340 (2002); Commun. Math. 
Phys. 229, 433 (2002). 

[10] S. Coles, An Introduction to Statistical Modeling of Extreme Values, Springer series in 
statistics (2001). 

[11] A. Soshnikov, Commun. Math. Phys. 207, 697 (1999); K. Johansson, Comm. Math. Phys. 

209, 437 (2000); I. M. Johnstone, Ann. Stat. 29, 295 (2001); P. Vivo, S. N. Majumdar, O. 

Bohigas, J. Phys. A: Math. Theor. 40, 4317 (2007) 
[12] A. Soshnikov, Electronic Communications in Probability 9, 82 (2004) and arxiv: 

|math.PR/0405090[ 
[13] G. Biroli, J. P. Bouchaud and M. Potters, EPL 78, 10001 (2007). 
[14] K. Johansson, Prob. Theo. and Rel. Fields 138, 75 (2007). 




(35) 



IV. CONCLUSION 



[15] M. Moshe, H. Neuberger, and B. Shapiro, Phys. Rev. Lett. 73, 1497 (1994). 
[16] O. Bohigas, J. X. de Carvalho, M. P. Pato, Phys. Rev. E 77, 011122 (2008). 
[17] O. Bohigas and M. P. Pato, Phys. Rev. E 74, 036212 (2006). 

[18] P. Forrester, Log-gases and Random Matrices www.ms.unimelb.edu.au to be published. 

[19] P. Forrester, Nonlinearity 19, 2989 (2006). 

[20] M. Dieng, Int. Math. Res. Notices 37, 2263 (2005). 

[21] A. C. Bertuola, O. Bohigas, and M. P. Pato, Phys. Rev. E. 70, 065102 (2004). 
[22] G. Akermann and P. Vivo, arxiv: math-ph/0806.1861vl. 
[23] G. Biroli, J. P. Bouchaud and M. Potters, Acta Phys. Pol. B 38, 4009 (2007). 
[24] O. Bohigas and M. P. Pato, Phys. Lett. B 595, 171 (2004). 
[25] S. P. Hastings and J. B. McLeod, Arch. Rat. Mech. Anal. 73, 31 (1980). 
Figure Captions 

Fig. 1 Density distribution of largest eigenvalue for GUE (J3 = 2). For complete 
sequence (/ = 1, Tracy- Widom) , for uncorrelated sequence (Weibull, limit / — > 0; Eq. 
( 135]) ) and partially incomplete sequences ( / = 0.3,0.6, from (|33|) ). 



